Supersymmetric Q-balls: A Numerical Study 
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We study numerically a class of non-topological solitons, the Q-balls, arising in supersymmetric 
extension of the Standard Model with low-energy, gauge-mediated symmetry breaking. Taking into 
account the exact form of the supersymmetric potential giving rise to Q-balls, we hnd that there is 
a lower limit on the value of the charge Q in order to make them classically stable: Q > 5 x 10 2 Q cr , 
where Q CI is constant depending on the parameters defining the potential and can be in the range 



1 < 



< 10 s 



If Q is the baryon number, stability with respect to the decay into protons 



requires Q > 10 17 Qcr, while if the gravitino mass is greater then m 3/ / 2 > 61MeV, no stable gauge- 
mediation supersymmetric Q-balls exist. Finally, we find that energy and radius of Q-balls can be 
parameterized as E ~ £eQ 3 ^ 4 and R ~ CkQ 1 ^ 4 , where £e and £r are slowly varying functions of 
the charge. 

PACS numbers: 05.45.Yv, 95.35.+d, 98.80.Cq 
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I. INTRODUCTION 

Cosmological and astrophysical observations indicate, 
almost undoubtedly, that our Universe is pervaded by 
a not-yet-known pressureless component, called "Dark 
Matter" [l| . Its energy density amount to about the 30% 
of the entire energy density of the universe. 
Promising non-baryonic, dark-matter candidates have 
been proposed since the discover of this mysterious com- 
ponent, such as axion or neutralino 

Recently enough ||, a particular class of non- 
topological solitons arising in supersymmetric extensions 
of the Standard Model, know as supersymmetric Q- 
balls Q, have been proposed as possible solution to the 
dark matter problem. 

In particular, Q-balls admitted in supersymmetric mod- 
els with low-energy, gauge-mediated symmetry break- 
ing |4|, are a plausible candidate for baryonic dark mat- 
ter |5|. 

Q-balls Q are lumps of matter, precisely a coherent 
state of a complex scalar field, carrying a conserved global 
charge. In the context of supersymmetric extensions of 
the Standard Model, the charge Q is some combination of 
baryon and lepton numbers, while the the scalar held is 
a gauge-singlet combination of squarks and sleptons cor- 
responding to some flat direction of the supersymmetric 
potential [3j]. In this class of models, supersymmetry is 
spontaneously broken at the scale Adsb ~ 10 7 GeV 
as a result, an effective potential for the flat directions 
arises 0| which, in turn, admits Q-balls as the non- 
perturbative ground state of the theory. 

A deep investigation ofgeneral properties [§] and as- 
trophysical implications [9( of Q-balls has been carried 
out in the last two decades. 

Up to now, however, the main properties of gauge- 
mediation supersymmetric Q-balls have been analyzed 



using an approximate expression of the potential giving 
rise to Q-balls. 

The aim of this paper is to study such a kind of Q-balls 
taking into account the exact form of the supersymmet- 
ric scalar potential quoted in Ref. 0]. We find that the 
expressions for energy and radius of Q-balls, which fully 
characterize their properties from a cosmological and as- 
trophysical viewpoint, can differ from the approximate 
case of about an order of magnitude. 



II. Q-BALLS: GENERAL PROPERTIES 

In this Section, we briefly review the Q-ball solution 
of a scalar theory with a global U(l) symmetry |6|. We 
consider a charged scalar field 4> whose lagrangian density 
is given by 



£=(g>^*)(c>V)-[/(M 2 



(1) 



In the next Section, we will identify 4> with one of the flat 
directions in supersymmetric extensions of the standard 
model, and will specify the form of the potential J7(|y|) 
which is not relevant in the present discussion. For the 
moment, we simply require its invariance under a global 
C/(l) symmetry. The corresponding conserved Noether 
charge q is normalized as 



q= i J (fix (ip* ip 



tptp 



(2) 



where a dot indicates a derivative with respect to time. 
(Throughout this paper, we follow the conventions of 
Ref. [lOj). For a given field configuration tp(t,r), the 
total energy is given by 



E 



d 6 x 



^| 2 + ^|Vyf 



(3) 
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We are interested to solutions of the field equations that 
correspond to a fixed value of the charge, namely Q, in 
Eq. ([2]). This is properly achieved by the introduction 
of the Lagrange multiplier lo associated to q, and by the 
requirement that the physical configuration makes the 
functional 



£uj — E 



(4) 



stationary with respect to independent variations of tp 
and ui |10| | . The requirement of time-independence of the 
total energy E implies the choice [1, Hoj 



<p(t,r) 



Atut 



<Kr), 



(5) 



with <p(r) real. Consequently, the functional £^ reduces 
to 



d 6 x 



1 



+ ojQ. (6) 
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The physical solutions have to satisfy the constraints 
6£u, „ 6£ 



= 0, 



5 cu 



= 0. 



(7) 



The first constraint leads to the equation of motion of 
the field Sir), 



dr 2 



dr 



dU 



(8) 



where, for simplicity, we assumed isotropy: 4> = </>(r) 
with r = \r\. 

The second constraint is equivalent to the requirement 
that the charge corresponding to the solution of the equa- 
tion of motion is equal to Q. 

A Q-ball is defined as the solution <p(r) of Eq. ||5J) satis- 
fying, at fixed charge Q, the boundary conditions <f>(r —> 
oo) = and d(j>/dr{r = 0) = |. 

In the next Section, we will analyze Q-ball configura- 
tions arising in a supersymmetric model where supersym- 
metry is broken via low-energy gauge mediation JiJ. 



III. SUPERSYMMETRIC Q-BALLS 

We consider a supersymmetric model in which super- 
symmetry is broken by low-energy gauge mediation 
In this kind of model the coupling of the massive vector- 
like messenger fields to the gauge multiplets, with cou- 
pling constant g ~ 10~ 2 , leads to the breaking of super- 
symmetry The coupling itself gives rise to an effec- 
tive potential for the flat direction tp whose lowest order 
(two- loop) contribution has been calculated in Ref. Q: 



U(z) 



A dx 
Jo 



z 2 — x{\ — x) + x(l — x) ln[x(l — x)z 2 



x(l~x)} 2 



(9) 
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FIG. 1: The potential U(<f>). Dotted lines refer to the asymp- 
totic expansions for small (<j> <C M) and large values (0 3> M) 
of the field <j> [see Eq. (fTD) ]. 



Here, z = <j>/M and M = M s /(2g), with Ms the mes- 
senger mass scale. The value of the mass parameter A 1 / 4 
is constrained as (see, e.g., Ref. p}): 10 3 GeV < A 1 / 4 < 
(g 1 / 2 /47r) v /TO 3/2 M P1 , where M P1 ~ 2.4 x 10 18 GeV is the 
reduced Planck mass and m%/2, the gravitino mass, is in 
the range 100 keV < m 3/2 < lGeV [EEI- 

The asymptotic expressions of U(z), for small and large 
z are Q: 



A 



(lnz 2 ) 2 -21nz 2 + 



if z < 1, 
if z > 1, 



(10) 



where m = V2A/M is the soft breaking mass and is of 
order ITeV p]]. In Fig. 1, we plot the potential U(z) 
with its asymptotic expansions (|10|) . 

Defining the critical charge Q cr = A/m 4 (whose mean- 
ing will be clear in the following), the constraint on A 
can be translated to a constraint on Q CI , namely: 



TeVV 



rn 



< Qcr < 10 8 



lio-v 



TeV V / m 3/2 
m ) VIOOkeV/ 



(11) 

A widely used approximation consists in replacing the 
full potential U{z) with its asymptotic expansions (JTUJ) 
in which a plateau plays the role of the logarithmic rise 
for large values of z. 

Within this approximation, it has been shown that the 
potential U(z) allows Q-balls solutions as the non per- 
turbative ground state of the model 0, [13] • Such states 
are known as supersymmetric Q-balls. In particular, 
for large charges, Q 3> Q cr , one can deduce analyti- 
cally the most important characteristics of Q-ball soli- 
tons. In more detail, the Q-ball profile is given by (3I. fT3| : 
<p(r) ~ 0o sm(ujr) / (tor) for r < R, and zero for r > R, 
where R = ir/u> is the radius of the Q-ball, and <fio = (j>(0). 
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FIG. 2: The Q-ball energy E normalized to the energy of Q 
quanta of the field <j> (that constitute the perturbative spec- 
trum of the theory) as a function of the charge. 



Moreover, one has [H, [13] 



Q 



m \Qc 
E 4V2tt / Q ^ 3/4 



-1/4 



TlQ cr 3 \Q, 

_R i_f n \V4 



V2 \Qc 



1/4 



(12) 
(13) 
(14) 
(15) 



Although the previous (simplified) analysis reveals the 
major properties of supersymmetric Q-balls, which are 
widely used in the literature in a cosmological and as- 
trophysical context, we wish to study them by taking 
into account the full potential ©. Since the form of the 
potential is involved, we need to solve the problem nu- 
merically. 

The computational procedure has been depicted in the 
previous Section. From a numerical viewpoint, it is sim- 
pler to fix uj rather than the total charge q. Then, once 
uj is fixed, we solve the equation of motion ([5]) with the 
condition d(f>/dr(r = 0) = 0. We look for the value of <f) 
such that the Q-ball solution exists; once this is achieved 
we insert 0(r) in Eqs. ^ and ([3]), obtaining the values 
of the charge and energy. Finally, we define the "radius" 
of the Q-ball, R, such that <p(R)/<f> = 0.1. 

In Fig. 2, we plot the quantity E/mQ as a function of 
the charge. If the energy E of the Q-ball at fixed charge 
Q is less then mQ, the soliton decays into Q quanta of 
the field (the perturbative spectrum of the theory) , each 
of them with mass m. Instead, if E < mQ the Q-ball 
is said to be classically stable, and then represents the 
ground state of the theory. 

Numerically, we find classical stability, E/mQ < 1, for 
Q > <3min, with Q min ~ 504Q cr . 

In Fig. 3, we plot w, the energy, the radius R, and </>o, 
as a function of the charge in the interval Q S [Q m in, 7.2x 
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FIG. 3: From top to down: the parameter ui, the energy E, the 
radius R, and (j>o as a function of the charge Q normalized to 
the critical charge Q CI . Dotted lines refer to the approximate 
solutions (fT6)) - (fT9)) . respectively [see also Eq. (f20)l and Table 
I]- 



l(F 7 Qcr]. 

If Q is the baryon number, stability with respect to 
the decay into protons requires the energy of a Q-ball 
to be less then E < m p Q, where m p ~ lGeV is the 
proton mass. We find numerically that this is attained 
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TABLE I: Nonlinear fit of the functions £'s defined in 
Eqs. (|16p - (|19p using a power-function of the type £(a;) = 
(a + bx p ) q . The parameters ei and £2 represent the max- 
imum percentage error of the functions £'s with respect to 
their numerical values in the 7.2 x 10 3r Qcr] 

and Q € [10 17 Q C r, 7.2 x 10 37 Q cr ], respectively. 





a 


b 


V 


q 


£1 


£2 


OJ 


- 9.567 


9.686 


0.381 


1 


1.32% 


0.16% 


E 


- 17.438 


15.559 


0.352 


1 


2.08% 


0.16% 


R 


- 7.162 


4.300 


0.560 


-0.712 


2.76% 


0.24% 


<j>0 


- 2.324 


1.292 


0.616 


0.668 


2.60% 


0.23% 



for charges larger then Q > 4.1 x 10 17 Q CI , where we 
assumed m = ITeV. 

We now wish to compare our numerical results to those 
obtained in the flat-potential approximation. Inspired by 
Eqs. (fT^ - tfTij)) . we write the quantities characterizing the 
Q-ball solution in the following way: 



,n = ^ ( l0gl ° 



Qc 



-1/4 



E 



£.E log! 



Q 



Q 



R 



1 = Zr[ log 



Qcr J \ Qc 

Q 



Qcr J V Qcr 



3/4 



1/4 



^0 

M 



u 1°, 



Q_ 

Qci 



1/4 



(16) 
(17) 
(18) 
(19) 



The functions £'s (which depend only logarithmically on 
the charge Q) parameterize the deviation from the simple 
power-laws (TT^1) - (|15I) . and are shown in Fig. 4. We fit 
their numerical values by the power-function 



£(x) = (a + bx p ) q . 



(20) 



In Table I, we report the values of the coefficients a, b, 
p, and q found by least-squaring the numerical data. We 
also show the maximum percentage error of the functions 
£'s with respect to their numerical values. In particular, 
£\ and £2 refer to the maximum percentage errors in the 
ranges Q € [Q ml „, 7.2 x 10 37 Q cr ] and Q G [10 17 Q cr , 7.2 x 
10 37 Q cr ], respectively. 

In is worth noting that, in the flat-potential approxi- 
mation, the functions £'s are constants whose values dif- 
fer from the numerical results of about an order of mag- 
nitude in the limit of large charges (say Q 3> 10 17 Q cr ). 

In Fig. 5, we show the Q-ball profile for three different 
values of the charge. We observe that to the charges 
Q = Qmin — 504Q cr and Q ~ 10 6 Q cr there corresponds 
the same value of radius, namely R ~ 11m -1 (see also the 
third panel in Fig. 3). However, looking at the shapes of 
the corresponding profiles, we see that in the first case the 



35 

30 
25 

2 20 

^3 

Vji IB 

10 
5 



10 15 20 25 30 35 
X 



X 

IS 20 



15 20 25 30 35 



10 15 20 25 30 35 
X 



15 20 25 30 35 
X 



FIG. 4: From top to down: the functions f w , £e, £r, and £4, 
defined in Eqs. (|16l) -(119 [ ). Dotted lines refer to the approxi- 
mating power-functions (|20[) (see also Table I). 



profile is spreader than the second one, indicating a larger 
"wall-thickness" of the Q-ball. Indeed, the larger the 
charge, the smaller is the wall-thickness, so that, as the 
charge increases the Q-ball approaches to the so-called 
"thin- wall" regime 0,0] (see continuous line in Fig. 5). 

As pointed out in Ref . 0] , for large values of the field 
4>, supergravity effects become important and give a con- 
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FIG. 5: Q-ball profile for three different values of the charge. 
Dotted line: Q — Q m m (R — 11m -1 ); Dashed line Q ~ 
10 6 Qcr (R ^ 11m -1 ); Continuous line Q ~ 7.2 x 10 37 g cr 
(i? ~ 3 x 10 8 m -1 ). 
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FIG. 6: The maximum charge Q max as a function of the grav- 
itino mass. The dotted line refers to the approximating func- 
tion (pTj) . 



tribution to the scalar potential of the form U gTav it y (4>) — 
m 3/2^ 2 ' approximatively. When this contribution domi- 
nates, E/gravity(0) 3> U (<fr) , Q-ball properties change dras- 
tically. Indeed, a different type of stable Q-balls are 
generated, the so-called "New- type Q-balls" The 
energy-charge relation is, in this case, -Enow typo ~ , m 3 / 2 Q 
(if Q is the baryon number, being m p > m 3 / 2 , the Q-ball 
is also stable with respect to the decay into protons). 
However, it is beyond the aim of this paper to study the 
numerical properties of such kind of Q-balls, and then 
we demand that the gauge-mediation potential §§§ dom- 
inates over the gravity-mediation one. 
Defining cq such that U'((f> cq ) = U' graVlty (</> oq ), it results 
U'{4>) > ^gravity (<A) f° r $ — 0eqi where a prime indicate 
differentiation with respect to <j). Here, following Ref. Q , 
we have compared the derivatives of gauge- and gravity- 
mediation potentials in order to determine their relative 
importance, since these are the quantity entering into the 
equation of motion ||5J). 

It is useful to introduce the "maximum charge" Q max 
such that, if Q < Q max then 4>o < </> C q- Since <j>Q = max0 
for a Q-ball configuration, the gravity effects can be ne- 
glected when Q < Q max . In Fig. 6, we plot <5 max as a 
function of the gravitino mass. It is easy to see that, 
neglecting logarithmic terms in the gauge-mediation po- 
tential, one finds U(<p) > [/gravity (0) for Q < Q max ~ 
(m/m 3 /2) 4 [where we used Eq. (|15p ]. Therefore, it is 
convenient to write the maximum charge as 

Qmax , f, m \ ( m \ 4 

= logio , (21) 

where £q(x) takes on the same form as in Eq. ([20)1 . By 
least-squaring the numerical data, we find a ~ 3.260, 
b ~ 5.750, p ~ —1.311, and q ~ 1.209, with a maximum 
percentage error on Q ma x of about 0.06%. Moreover, 
we find Q m ax ^ 4.9 x 10 28 and Q max ^ 6.4 x 10 12 for 
m 3/2 — 100 keV and m 3 / 2 = lGeV, respectively. 

If Q is the baryon number, it is interesting to observe 



that, depending on the value of the gravitino mass, it can 
be Q max < 4.1 x 10 17 Q cr , indicating that no stable Q-ball 
solution exists. Numerically, we find the this happens for 
gravitino masses above m 3 / 2 ^ 60.8MeV. 

IV. CONCLUSIONS 

In this paper, we have studied Q-balls-type solutions 
admitted in supersymmetric particle-physics models with 
low-energy, gauge-mediated supersymmetry breaking. 
Taking into account the exact form of the supersymmet- 
ric potential, we have analyzed classical stability of Q- 
balls. We have found, numerically, that only Q-balls with 
charge Q > 5 x 10 2 Q cr are stable against the decay into 
quanta constituting the perturbative spectrum of the the- 
ory. Here, the "critical charge" Q C i is a model-dependent 
parameter given in Eq. (|11|) . 

Moreover, if the conserved charge Q is the baryon num- 
ber, stability with respect to the decay into protons (the 
lightest baryonic particle) requires Q > 10 17 Q C r- 

Although no analytical expressions for the quantity 
characterizing Q-ball solutions can be found, we were 
able to approximate the numerical results by suitable 
functions: we have found, indeed, that energy and ra- 
dius of Q-balls, which fully characterize their properties 
from a cosmological and astrophysical viewpoint, can be 
parameterized as E ~ £,eQ 3 ^ 4 and R ~ ^rQ 1 ^ 4 , where 
(,e and £,r are slowly varying functions of the charge [see 
Eqs. (fT7 ]) -(|T8| l . Eq. ([20]), and Table I]. 
In the (approximate) case of exactly flat potential con- 
sidered in the literature, the functions £'s are constants 
whose values can differ from our results of about an order 
of magnitude. 

For large values of the scalar condensate defining a 
supersymmetric Q-ball, supergravity effects become im- 
portant and give a contribution to the scalar potential 
which, in turn, change drastically the Q-ball properties. 
Introducing the "maximum charge" Q max such that, if 
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Q < Qmax supergravity effects are negligible, we have 
found that <5 max , as a function of the gravitino mass, 
can be well approximated by Q max ~ £q to 3~/2> wl ^ a £q a 
slowly varying functions of m 3 / 2 [see Eq. (|2Tj)]. 
In particular, if Q is the baryon number we have found 
that, for gravitino masses above m 3 / 2 ^ 61MeV, the 
maximum charge is less then the charge required for clas- 
sical stability, indicating that no stable gauge-mediation 
supersymmetric Q-balls exist. 
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